Kinetics of condensation of flexible polyelectrolytes in poor solvents: effects of solvent 

quality, valence and size of counterions 
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The collapse kinetics of strongly charged polyelectrolytes in poor solvents is investigated by Langevin 
simulations and scaling arguments. The rate of collapse increases sharply as the valence of counte- 
rions, z, increases from one to four. The combined system of the collapsed chain and the condensed 
counterions forms a Wigner crystal when the solvent quality is not too poor provided z > 2. For very 
poor solvents the morphology of the collapsed structure resembles a Wigner glass. For a fixed z and 
quality of the solvent the efficiency of collapse decreases dramatically as the size of the counterion 
increases. A valence dependent diagram of states in poor solvents is derived. 



PACS numbers: 36.20.-r, 61.41+e 

Upon condensation of DNA into toroidal structures its 
volume decreases by nearly four orders of magnitude jlj . 
Counterion condensation, which leads to attractive inter- 
actions between the segments of the polyanion chain, is 
a key factor in determining DNA collapse jl. Experi- 
ments show that counterions with trivalent (or greater) 
valence, z, are often required to induce chain collapse 
0. Precipitation (requiring chain contraction) in aque- 
ous solution of a highly charged flexible polyelectrolyte 
(PE) (polystyrene sulfonate) occurs only when counte- 
rions with z = 3 or 4 are added B. The collapse of 
DNA (a stiff chain) depends on the interplay of several 
effects such as bending rigidity, valence and size of coun- 
terions, and hydration forces Since intrinsic stiff- 
ness is not relevant for flexible PE the collapse transition 
may be simpler to describe than in stiff PE. However, 
because of competing effects of the bare electrostatic in- 
teractions, counterion condensation and the solvent qual- 
ity the phase diagram of highly charged flexible PEs is 
quite complicated 0]. Simulations |3]||, scaling argu- 
ments 0,0,^1, and theory || are beginning to provide a 
picture of the possible structures that emerge in highly 
charged polyelectrolytes in which correlations induced by 
counterions play a crucial role fl-|r^,p|,|io| . 

In this letter we describe the dynamics of collapse of 
strongly charged flexible PE in poor solvents with em- 
phasis on the mechanisms of collapse and their depen- 
dence on the valence, size of counterions and solvent qual- 
ity. We consider polyelectrolyte chains that are highly 
charged so that f(ls/b) 2 > 1 where f(=l in the exam- 
ples considered here) is the fraction of charged monomers, 
Ib = e 2 /Ai:ekBT is the Bjerrum length, with T being 
temperature, e is the dielectric constant of the solvent, 
and b is the size of the monomers. The major results 
of this study, which are obtained using a combination of 
Langevin simulations and scaling arguments, are: 

(a) The mechanisms and kinetics of collapse dramat- 
ically depend on the valence z of the counterions. For 
all range of parameters the slowest kinetics is found for 
z = 1, and generally the rate of globule formation in- 



creases with z. 

(b) The equilibrium structure of the collapsed poly- 
electrolyte depends critically on the solvent quality. If 
the solvent is not too poor then the conformation of the 
counterion and the charged polymer is a Wigner crystal. 
For z > 2, we find a BCC-like crystal. When the sol- 
vent quality is very poor amorphous structures (Wigner 
glasses) are formed. 

(c) For a fixed solvent quality, the efficiency of chain 
collapse depends dramatically on the size of the counte- 
rions. The rate of approach to the globular conformation 
increases as the size of the counterion decreases at a fixed 
value of z. 

The polyelectrolyte is taken to be a flexible chain con- 
sisting of N monomers each with a charge of —1 (mea- 
sured in units of e). Successive monomers are connected 
by an elastic spring with spring constant 3ksT/b 2 . The 
valence of the cations varies from z = 1 — 4 and their 
density is determined by the condition of charge neu- 
trality. The non-electrostatic interaction between the 
particles (monomer or counterion) i and j with radii ai 
and <jj respectively that are separated by ry is taken 
to be H LJ ( riJ )/k B T = £lj[(^) 12 - 2{^-f] where av- 
erage bond length r Q = a% + <7j . The quality of the sol- 
vent is expressed in terms of the second virial coefficient 

H LJ 

v 2 = L d 3 r(l — e ^bt) which is negative for poor sol- 
vents. 

The chain and the counterions are placed in a trun- 
cated octahedron (Wigner— Seitz cell of BCC) unit cell 
with the lattice constant L and the volume of the unit 
cell is V = T^L 3 . The long range nature of the Coulomb 
interaction is treated by accounting for interactions with 
the image charges in a periodically replicated system p| . 
The computations were done with N = 120 (L — 40cr m ) 
or N = 240 (L = 80er m ). Most of the results presented 
here are for the latter case so that the monomer density 
p m = N/V = 9.36 x 10- 4 cr m 3 , where a m (= 6/2) is the ra- 
dius of the monomer. We start from either 0-conditions 
or from high temperature and quench the temperature 
below Manning JL2| condensation condition. We have 
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verified that the results (collapse mechanisms and the 
morphology of the collapsed structures) do not depend on 
the initial conditions. The chain dynamics is described 
by the Langevin equation. The equations of motion are 
integrated with a step size = 10~ 3 r where the unit of 
time r = b 2 /D — b 2 (/ksT, where £ is the friction co- 
efficient. In order to obtain results with relatively small 
errors averages of various quantities are taken over at 
least 40 independent initial conditions. 
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FIG. 1. The time dependence of the radius of gyration 
following a quench from 0-solvent to poor solvent condi- 
tion. The value of N = 240 and l B = 5.36. Fig. (la) is for 
v 2 = -0.066 3 where as Fig. (lb) is for v 2 = -3.696 3 . The 
numbers on the curves denote the valence of the counterions. 

The time dependence of < R g (t) > / R c for various 
values of z, where R c is the size of the compact globule, 
shows that for v 2 = -0.066 3 ,Z B = 5.36 (Fig.(la)) only 
counterions with z = 3 or 4 induce globule formation 
on the time scale of the computations. The structure of 
the chain with z = 2 is extended at t ~ 300r. When 
the quality of the solvent is poorer (v-2 = — 3.696 3 ) then 
the condensation occurs as long as z > 2 (see Fig. (lb)). 
Monovalent counterions are inefficient in causing chain 
collapse. The rate of approach to globular conformation 
is largest for z = 4 and it decreases as the valence is 
reduced. Although the final equilibrium structure is ex- 
pected to be the same, the kinetics is sensitively depen- 
dent on the valence of the counterions. For z > 3, the 
difference in kinetics becomes negligible as the solvent 
quality decreases. 

The approach to the globular structure occurs in 
roughly three stages for all z. In the initial stage coun- 
terions condense onto the polyanion. This takes place 
on a diffusion controlled time scale r cc ss pm^/D which 
for our simulation is « 25t. We calculated the time re- 
quired for 90% of counterions to condense. The time 
scales are determined only by p m and are independent of 



z and v%. In the second stage globular clusters contain- 
ing both monomers and counterions form. The clusters, 
the number of which is dependent on z and are con- 
nected by "strings" giving rise to pearl-necklace struc- 
tures (see below). This process occurs in the time regime 
t cc < t < tclust- For t > tclust, coarsening of the 
clusters takes place till the equilibrium globular structure 
is formed jljj. 

The variations in the collapse kinetics with z may be 
understood in terms of formation of clusters (or pearls) 
and subsequent growth (or merging) of such clusters into 
a compact globule H]. When monovalent counterions 
are condensed they combine with the charges on the 
backbone of the polyanion to form ion pairs. Even though 
the condensed counterions are mobile along the polyanion 
backbone, upon condensation randomly oriented dipolcs 
of magnitude p w eb are formed. In order to form con- 
tacts between the segments of the chain the attractive 
energy between the dipoles must exceed &bT. Thus, 
the range of attractive interactions r t w (/^fr 2 ) 1 / 3 < Ib- 
Since rt/ls < 1 we expect cluster formation to occur 
only on short length scales. This generates a large num- 
ber of clusters and their subsequent merger (nucleation) 
to form a fully compact globule is expected to be slow. 
In fact, we do not see globule formation with z = 1 on 
the time scales of our simulation for any value of V2- 

A very different mechanism emerges when chain col- 
lapse is induced by multivalent cations. When the Man- 
ning parameter = ls/b > 1/z some of the ions 
are condensed [0]. A counterion can locally neutral- 
ize the negative charge on the monomer, and the re- 
sulting bound species has an effective charge (z — l)e. 
Bare monomer charges that are separated by a large dis- 
tance along the contour can be attracted to the positive 
charge - a process referred to as "ion-bridging" |3|]. The 
range over which such attractive interactions are effec- 
tive clearly increases with z. Since the size (and hence 
the number of clusters ) of the clusters is controlled by 
the length of the contour over which ionbridges form, 
it follows that the efficiency of collapse should increase 
with z. The time dependence of the correlation function 

c(t) = Ei <j ,\i- j \> 5 Q ( 2b - \n - r j\) > where O) is thc 
step function, shows that for z = 2, 3 and z = 4 con- 
tacts between monomers that are separated by at least 
four bonds occur more rapidly as z increases. The ion- 
bridging mechanism leads to rapid condensation of the 
chain to metastable compact pearl-necklace structures for 
z = 3 and 4 Jl|@. 

For z > 2 the equilibrium globular conformation is 
kinetically reached provided that the solvent is suffi- 
ciently poor. However, the arrangement of the monomers 
and counterions in the globular structure depends on 
V2- From Fig. (2a), which shows a plot of the time 
development of the pair distribution function g(r) for 
i>2 = — 0.346 3 , Ib = 66 and z = 2, we infer that the com- 
bined system of counterions and the polyelectrolyte forms 
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a Wigner crystal jlTj . The formation of Wigner crystal is 
seen by the presence of several peaks in g(r) at long times 
(t > 350t). The analysis of the peak position reveals that 
the ordered structure corresponds to a BCC lattice. The 
ordering is considerably fuzzy (see Fig. (2b)) when the 
solvent quality is extremely poor (i>2 = — 62.346 3 ). The 
resulting structure is amorphous (Wigner glass). Thus, 
for «2 = — 62.346 3 the dynamics of approach to the or- 
dered state is slow i.e., there are manifestations of glassy 
dynamics. 
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FIG. 2. Pair correlation function for the combined system 
of the PE chain (N — 240) and the counterions (z = 2) in long 
time limit (t > 350r). The system is quenched to Ib = 66. 
The top panel is for «2 = — 0.346 3 where the bottom panel 
corresponds to V2 = —62.34b 3 . The insets show the corre- 
lation functions between monomers (symbols) and between 
counterions (solid lines). The lack of correlations between 
counterions is evident. 

These observations can be used to construct a valence 
dependent phase diagram for strongly charged polyelec- 
trolyte chain in poor solvents Jl8| . Upon counterion con- 
densation the total charge of the polyelectrolyte decreases 
from Nf to Nf « k(L/Ib)(1/z) where k w — ln<f) and <j> 
is the volume fraction of the free counterions 0. The 
size of the polyelectrolyte is given by L w n 2 b 2 N /IbZ 2 
provided Ib > K^z^bN 1 ' 2 . As the quality of the sol- 
vent decreases to a level such that thermal blob size 
£t ~ k 4 /l w 2| < £e/ « Ibz 2 /k 2 (size of the electrostatic 
blob) then the chain condenses to a globule. The bound- 
ary dividing the stretched and collapsed conformation is 
given by |t>2| ~ b 4 K 2 /Ibz 2 . In the globular phase we find 
two regimes, one corresponding to the Wigner crystal 
and the other a Wigner glass. The boundary between 
the two is obtained by equating the gain in the energy 
upon condensation (f» ze 2 /de) to the attractive interac- 
tion due to the poor solvent quality (kTv 2 b~ 6 ). This 
leads to the condition jt^l ~ z l / 2 b 5 / 2 lg. To validate 



the phase diagram we performed extensive simulations 
for all z(= 1,2,3,4) at valuase of Ig 1 and \v 2 \ indicated 
by asterisks in Fig. (|J) . The structural conformations are 
in accord with the predicted phase diagram. 
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FIG. 3. Valence dependent diagram of states in the (|i)2| 
and Ig 1 ) plane for strongly charged PE. The dashed lines rep- 
resent the boundary between stretched and collapsed states 
and depend on z. The z-dependent solid lines in the col- 
lapsed region separate Wigner crystalline region from the 
Wigner "glassy" region. The asteriks are the simulation re- 
sults with each data point corresponding to four z values 
(z = 1—4). Pictures of the conformation of the chain 
in region A,B and C are also shown. The lg values are 
indicated by arrows. The letter (a-1) near the asterisks 
corresponds to values of \v 2 \/b s = 62.34,62.34,25.82,25.82, 
15.20,12.36,7.51,4.47,3.91,3.69,0.34,0.06 respectively. 

Oosawa has argued that the fraction of condensed 
counterions is not significantly changed by altering the 
size of the counterions provided the volume fraction of 
the polyanion is small . Since the correlation between 
the monomers and counterions is dependent on the size of 
the counterions the dynamics of collapse should be sensi- 
tive to R = <y m j '<t c where a c is the size of the counterions. 
We carried out simulations for R = 2 and R = 5/7 with 
v 2 = -3.696 3 for N = 240 at l B = 5.36. For R = 2 we 
find that the rate of globular formation in enhanced com- 
pared to R = 1 for z — 2,3 and 4. The difference between 
the two cases is most dramatic for z = 2 (see Fig.(^)). For 
the divalent cation globular formation occurs at t f=s 300t 
when R = 1 (see Fig. (la)). However, when the size of 
counterions is made smaller the PE chain almost reaches 
the globular state at t w 150r. Conversely, when the size 
of the counterion is increased the efficiency of collapse 
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decreases (see Fig.(Q)). This establishes the importance 
of the size of the counterfoil, and in particular the corre- 
lation between them in determining collapse dynamics. 
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FIG. 4. The time dependence of < R g (t) > (N = 240) 
for different sizes of the counterions (z = 2) for 
V2 — —3.69b 3 ,1b = 5.36. The top curve (dotted line) is for 
R = 5/7 and the dashed line is for R — 2. For comparison 
results for the condition R — 1 are also displayed as solid line. 

The findings presented here show that a bewildering 
range of dynamical behavior is to be expected in the col- 
lapse of strongly charged polyelectrolytes mediated by 
countcrion condensation. The prediction that the for- 
mation of ordered states Wigner crystals and amorphous 
structure (Wigner glasses) should depend on the quality 
of the solvent should be amenable to experimental (neu- 
tron or light scattering) tests. Even if the Wigner glass is 
metastable our theory suggests that the dynamics in the 
two regions (see Fig.(||)) covering the globular states of 
PE should be dramatically different and depend on the 
valence of the counterions. The prediction that size con- 
trols the kinetics of collapse can explains the changes in 
the efficiency of DNA condensation with differing cations 
with fixed valence |ll| . 
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